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EFFECTS  OF  OCEAN  FRONTS  ON  MULTIPATH  VERTICAL  ANGLES 
AND  THE  ASSOCIATED  BIASES  IN  ARRAY  BEARING  ESTIMATION 

R.P.  Flanagan  and  X.  Zaba lgogeazcoa 


An  oceanic  frontal  system  perpendicular  to  tne  propaga¬ 
tion  path  is  modeled  with  a  deep  receiver  in  warm  water  and 
a  shallow  source  which  is  moved  through  the  frontal  system 
into  progressively  colder  water.  A  ray  model  is  used  to 
calculate  the  arrival  structure  and  determine  the  effects 
of  the  frontal  system  on  ray  amplitude  and  vertical  angle. 


The  receiving  horizontal  line  array  determines  the  bias  in 
bearing  estimation  resulting  from  the  changing  multipath 
structure.  This  bias  is  found  to  be  heavily  dependent  on 
array  orientation  and  range  into  the  frontal  system  while 
being  only  slightly  dependent  on  array  beamwidth.  The  under¬ 
lying  physics  causing  the  multipath  variations  is  the  frontal 


system  acting  as  a  focusing  lens.  When  the  transducer  posi- 


ions  arc  reversed,  that  is  a  deep  receiver  in  cold  water 


Introduction 


The  propagation  of  sound  through  oceanic  fronts  is 
a  topic  of  interest  due  to  the  frequent  appearance  of 
these  systems  in  the  world  oceans.  Typically  an  ocean 
front  may  be  visualized  as  a  well  defined  transition 
region  between  two  volumes  of  water  with  characteristic 
sound  velocity  distributions.  One  of  the  effects  of 
such  a  region  in  a  propagation  path  is  the  alteration  of 
the  multipath  structure  of  the  received  signal;  in  fact 
an  oceanic  front  may  act  as  an  acoustic  lens  with  either 
convergent  or  dispersive  characteristics  depending  on  the 
direction  of  transmission.  Since  source  bearing  estima¬ 
tion  with  a  linear  array  is  highly  dependent  on  the 
vertical  structure  of  the  acoustic  paths  at  the  sensors 
one  of  the  major  effects  of  an  oceanic  front  is  the  in¬ 
troduction  of  biased  bearing  estimation.  The  object  of 
this  paper  is  the  determination  of  the  magnitude  of  this 
bias  for  a  typical  frontal  system  anomaly. 

In  this  study  ray  tracing  techniques  applied  to  an 
ocean  with  vertical  and  horizontal  variations  in  the  sound 
velocity  profile  have  been  used  to  analyze  the  effect  of 
a  frontal  system  on  acoustic  transmission.  The  bulk  of 
the  material  is  concerned  with  a  source  located  in  pro¬ 
gressively  colder  water  propagating  a  200  Hz  CW  signal 
through  the  frontal  system  to  a  warm  water  array.  After 
the  multipath  structure,  arrival  angle,  bearing  estimation 
and  bias  are  thoroughly  presented  for  this  transducer 
geometry  the  positions  of  the  source  and  receiver  are 
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reversed  and  a  brief  examination  of  the  frontal  effects 
is  made.  In  both  cases  the  bias  in  bearing  estimation 
introduced  by  the  frontal  system  was  found  to  be  a 
function  of  ranqe,  beamwidth,  and  directional  orienta¬ 
tion  of  the  array. 
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I  *  The  Frontal  System 

An  oceanic  front  may  be  defined  as  the  transition 
region  which  divides  water  masses  of  different  physical 
characteristics  with  the  surface  positions  of  the  fronts 
being  marked  by  stiong  horizontal  temperature  gradients. 
These  fronts  tend  to  remain  in  one  area,  but  can  exhibit 
severe  variations  within  that  area.1-2  A  general  outline 
of  the  frontal  systems  of  the  world  oceans  along  with 
their  seasonal  positions  is  given  by  Laevastu  and  LaFond.^ 
Work  on  specific  systems  includes  Levenson  and  Doblar*  for 
the  North  Atlantic,  Roden^  for  the  North  Pacific  and 
Bannister,  et.  al.6  for  the  South  Tasman  Sea. 

The  system  we  have  chosen  to  model  is  bi-frontal 
and  extends  from  warm  to  cold  water  regions,  (Fig.  1). 
Proceeding  from  warm  to  cold  water  the  first  region  has 
a  profile  with  a  negative  gradient  extending  from  the 
surface  to  the  SOFAR  axis.  The  profile  for  the  second 
region  has  a  near  isovelocity  layer  extending  from  the 
surface  to  150  m.  Beyond  that  depth  the  sound  velocities 
gradually  merge  to  form  an  elevated  SOFAR  axis.  The  area 
separating  these  two  regions  is  termed  the  first  front 
and  is  50  km  wide.  The  third  region  contains  the  coldest 
water  and  has  a  near  isovelocity  layer  in  the  upper  200  ■ 
joined  to  a  positive  gradient  below  that  depth.  The 
second  front,  which  separates  regions  2  and  3,  is  also 
50  km  wide.  Because  of  the  imaginary  location  of  the 
receiving  array,  and  thus  its  variability  in  range,  it 
is  possible  to  examine  the  propagation  of  sound  for  both 
the  single  and  double  front  cases. 


The  top  half  of  Fig.  1  shows  the  location  of  the 
fronts  with  range.  The  receiver  is  at  2100  m  in  the 
warm  water  profile  while  the  100  m  deep  source  was 
placed  at  ranges  varying  from  700  to  1125  km.  Each 
region  has  a  range  independent  sound  velocity  while 
the  fronts  contain  the  transition  profiles  between 
surrounding  regions. 

The  fronts  are  taken  to  be  perpendicular  to  the  pro¬ 
pagation  path.  If  this  is  not  the  case  horizontal  refrac¬ 
tion  takes  place  and  the  rays  propagate  out  of  the  vertical 
plane  indicating  a  source  of  bearing  bias  which  will  not 
be  examined  here.  Finally,  the  width  of  the  front  will 
affect  the  bias  in  bearing  estimation.  Calculations  with 
a  front  width  of  5  km  showed  rapid  oscillations  of  this 
bias  with  range.  This  occurs  because  the  rays  alternately 
sample  and  miss  the  front  since  the  anomaly  decreases  in 
amplitude  with  depth  and  the  rays  may  therefore  pass  under 
it  at  certain  ranges. 


The  ray  model  used  for  this  study  is  described  else¬ 
where7  and  is  capable  of  tracing  rays  in  a  medium  with 
vertical  and  horizontal  variations  in  the  sound  speed  pro¬ 
file.  The  basic  profile  consisted  of  four  pieces  with  a 
linear  bottom  section,  two  hyperbolic  cosines  for  the  SOFAR 
region,  and  a  surface  section  with  a  sound  velocity  distribu¬ 
tion  that  yields  a  conic  ray  path.  The  horizontal  variation 
of  the  profiles  occurred  only  in  the  two  frontal  regions 
with  ranoe  steps  of  ICO  m.  The  bottom  was  assumed  to  be 
flat  and  lossy  to  limit  the  significant  arrivals  to  refracted- 
refracted  (RR)  and  ref racted-surf ace  reflected  (RSR)  rays. 

The  corresponding  profiles  and  profile  anomalies  are  shown 
in  the  lower  half  of  Fig.  1  and  Fig.  2  respectively.  It 
can  be  seen  that  the  anomalies  for  these  two  fronts  are 
similar  in  general  form. 


III.  The  Array  Reception  Model 

The  receiving  array  was  modeled  as  a  horizontal 
linear  set  of  phones.  In  order  to  avoid  specifying  the 
number,  separation  distances,  and  shading  of  the  phones, 
the  calculations  were  performed  for  3  dB  beamwidths  of 
1°,  2°  and  5°  with  sin  x/x  weightings.  The  range  of 
beamwidths  was  based  on  the  examination  of  several  papers 
on  arrays  a  few  of  which  are  ref.  8-11.  Typically  beam- 
width  broadens  when  passing  from  broadside  to  endfire 
aspects  so  the  results  for  the  different  beamwidths  can 
be  mixed  to  more  accurately  depict  this.  The  array  recep¬ 
tion  as  a  function  of  array  axis  orientation  with  respect 
to  the  frontal  system  was  simulated.  The  array  axis  was 
rotated  from  a  position  parallel  to  the  fronts  to  the 
perpendicular  position  in  1°  increments.  For  each  array 
orientation  the  array  response  as  a  function  of  azimuth 
(9)  was  obtained  by  means  of  the  incoherent  beamformer. 


where  A  and  k  arc  the  amplitudes  and  vector  wave  numbers 
n  n 

of  the  N  paths  determined  by  the  ray  tracing  model  ,yti  n 
is  the  vertical  elevation  of  the  paths,  kQ  is  the  horizontal 
array  steering  vector  and  L  is  a  parameter  that  adjusts 
the  desired  beamwidth. 

The  geometry  of  the  array  is  shown  in  Fig.  3.  The 
main  lobe  corresponding  to  a  steering  to  azimuth  0Q  is 
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shown  as  a  figure  of  revolution  about  the  array  axis. 

If  a  source  is  located  at  azimuth  0^  only  rays  reaching 
the  array  with  vertical  arrival  angles  between^^^and 
^f^will  be  processed  in  the  beamformer  output. 

Figures  4  and  5  show  vertical  angle  of  arrival  vs 
bearing  (0°  being  endfire)  for  1°  and  5°  beamwidths. 

From  these  plots  it  is  possible  to  obtain  a  bearing 
estimate.  The  two  horizontal  lines  on  Fig.  5  indicate 
arrivals  with  vertical  angles  of  7.5°  and  15°.  Let  the 
vertical  line  at  20°  bearing  represent  the  true  bearing 
of  the  source.  The  intersection  of  these  lines  at  points 
A  and  D  may  be  continued  to  the  bearing  axis  by  the 
dashed  lines.  Apparent  bearings  of  21.25°  and  25°  yield 
biases  of  1.2°  and  5°  for  the  two  rays.  If  equi- amplitude 
rays  are  chosen  then  an  averaging  process  can  be  used  to 
set  the  bias  at  3.1°.  If  the  true  bearing  is  50°  the 
points  C  and  D  convert  to  biases  of  .3°  and  1.8°  with 
an  average  of  1.0°.  From  the  comparison  of  these  two 
examples  it  can  be  seen  that  this  bias  is  a  function  of 
array  orientation  with  the  largest  biases  occurring  in 
the  endfire  aspect. 

If  ray  amplitudes  are  different  another  m-  thod  of 
determining  the  bias  in  bearing  estimate  is  needed. 

Those  characteristics  necessitated  the  incoherent  sum¬ 
mation  of  the  rays  over  a  specified  beamwidth  and 
weighting  function  described  previously.  A  three 
dimensional  plot  of  array  energy  vs  true- bearing  from 
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array  vs  apparent  bearing  from  array  yields  a  visual  repreaen- 
tacion  of  the  bias  in  bearing  estimation  as  will  be  shown  for 
the  specific  front  process. 


The  Computer  Simulation 


The  ray  model  was  run  at  various  ranges  before  the 
fronts  and  in  the  fronts.  The  ranges  within  the  frontal 
region  were  also  run  without  the  frontal  effects  in  order 
to  examine  the  changes  that  the  fronts  had  caused.  All 
of  the  data  was  examined  with  beamwidths  of  1°,  2°  and  5°. 

As  an  example,  three  dimensional  plots  and  the  re¬ 
sultant  bearing  bias  vs  true  bearing  diagrams  are  given 
for  ranges  of  875  and  910  km  and  a  1°  beamwidth.  Figure  1 
shows  these  two  ranges  to  be  in  the  center  of  and  10  km 
after  the  first  front  respectively.  Figures  6  and  7 
represent  the  array  response  for  these  ranges  without 
the  fronts,  and  Figs.  8  and  9  depict  the  complementary 
data  but  include  the  frontal  effects.  Figure  10  depicts 
the  bias  in  bearing  for  each  of  the  four  cases  above. 
Examining  Fig.  6  with  a  bearing  of  source  from  array  of 
0° ,  three  peaks  are  observed  at  apparent  bearings  of  9.5°, 
11.1°,  and  13.6°.  These  peaks  represent  the  energy  from 
ray  pairs,  each  pair  with  the  same  number  of  cycles 
between  the  source  and  receiver  and  consisting  of  a 
positive  and  negative  angle  at  the  source.  It  is  ap¬ 
parent  that  a  horizontal  array  in  the  endfire  aspect, 
with  a  beamwidth  and  directional  increment  of  .1°  can 
separate  ray  arrivals  with  vertical  angle  differences  of 
twice  the  beamwidth  if  the  array  operates  as  presumed. 

As  the  array  is  moved  towards  broadside  the  ability  to 
separate  ray  arrivals  is  lost.  This  is  observable  from 


Fig.  4  as  the  range  of  vertical  angle  increases  for 
constant  bearing  as  bearing  is  increased.  Returning  to 
Fig.  6  we  define  the  apparent  bearing  (bearing  estimate) 
of  the  source  frcm  the  array  as  the  point  of  maximum 
energy;  then  11.  1 4  is  the  bearing  bias.  Naturally  it 
is  desirable  to  make  the  bias  as  small  as  possible  and 
we  would  like  to  pick  the  peak  at  9.5  .  This  could  be 
acomplished  if  we  used  the  first  deviation  from  zero 
as  apparent  bearing  is  increased,  however  it  nay  not 
bo  possible  to  determine  if  it  is  single  or  multiple 
sources  creating  the  multi-bearing  readings. 

Comparing  Fig.  6  with  8  and  Fig.  7  with  9  wo  can 
examine  the  effects  of  the  first  front.  There  appears 
to  be  more  energy  for  the  frontal  cases,  and  it  is 
shifted  towards  the  lower  end  of  the  apparent  bearing 
axis  reducing  the  bearing  bias.  This  is  also  indicated 
in  Fig.  10.  Figures  11  and  12  show  the  amplitudes  of 
the  positive  source  angle  rays  which  reach  the  array 
with  vertical  angle  for  the  cases  without  the  front 
and  with  the  front  respectively.  The  arrival  sets  are 
shown  for  different  ranges  and  the  location  of  these 
points  with  respect  to  the  frontal  regions  can  be  deter¬ 
mined  from  Fiq.  1.  Fxamining  the  two  Figs,  at  875  km 
and  910  km  there  is  a  definite  increase  in  the  number 
of  arrivals  and  a  decrease  in  8^  for  the  arrivals  when 
the  effects  of  the  frontal  region  are  included. 
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One  simple  case  is  presented  in  the  Appendix  showing 
transmission  from  cold  to  warm  water  through  a  front  focuses 
the  multipath  structure.  That  is,  rays  from  the  cold  water 
which  interacted  with  the  surface  will  not  reach  the  surface 
in  warmer  waters.  Therefore,  the  acoustics  with  the  cold 
to  warm  water  front  and  these  approximate  source  and  re¬ 
ceiver  depths  should  support  more  RR  type  rays  with  higher 
cycle  number  and  lower  SR  than  the  acoustics  without  the 
f  ront . 

Figure  13  depicts  the  arrival  pattern  for  the  positive 
source  angles  before  the  fronts  from  750  km  to  850  km.  It 
is  possible  to  detect  trends  in  this  pattern.  If  the  bases 
of  the  rays  are  joined  in  the  760  km  to  810  km  ranges  a 
shape  is  formed  similar  to  the  greek  letter  u  with  the 
final  tail  clipped  off.  Figure  14  is  an  expansion  of  the 
previous  figure  including  ranges  from  750  km  to  1100  km 
without  the  frontal  process.  Each  u  shape  here  represents 
rays  with  the  same  cycle  number  as  indicated  on  the  figure. 
Because  there  are  no  horizontal  variations,  i.e.,  profile 
1  is  used  throughout,  the  angles  may  be  found  where  pro¬ 
pagation  type  changes  from  RR  to  RSR  to  SRBR.  These  divi¬ 
sions  are  indicated  on  the  plot.  A  rough  approximation  of 
spreading  loss  can  be  determined  by  examining  dR/d0R  of 
each  curve.  Large  values  of  this  derivative  indicate  low 
amplitudes  while  small  values  indicate  large  amplitudes. 
Therefore,  the  amplitudes  of  the  RR  ray  types  increase  as 
0R  decreases.  The  RSR  rays  have  low  amplitudes  at  the 
transition  regions  and  large  amplitudes  in  the  "valley" 


regions  indicated  in  the  figure.  Further,  because 
dR/dOR  is  larger  there,  the  RSR  rays  or.  the  left  side 
of  the  valley  have  lower  amplitude  than  those  on  the 
right.  The  SRBR  rays  have  consistently  high  losses  due 
to  spreading  which  are  further  magnified  by  the  bottom 
losses,  so  these  rays  will  not  be  mentioned  any  further. 
Using  this  figure  it  is  possible  to  determine  the  number 
and  type  of  ray  solutions  at  a  given  range.  For  example, 
at  900  km  the  arrival  set  consists  of  a  16  cycle  RR  ray 
and  two  16  cycle  RSR  rays.  At  1040  km  the  arrival  set 
contains  one  HR  and  two  RSR  rays  with  18  cycles  and  the 
same  types  with  19  cycles. 

The  bias  in  bearing  estimation  can  be  examined 
using  this  figure  and  the  general  lule  that  an  increase 
in  vertical  angle  causes  an  increase  in  bearing  bias. 
Starting  at  750  km  and  increasing  range  the  angle  of 
the  arrival  with  the  largest  .amplitude  at  each  ranqe 
will  indicate  the  bearing  bias.  From  750  km  to  850  km 
at  10  km  Intervals  th*>  angler  are:  8.8,  10.5,  8.2,  8.3, 
8.5,  9.7,  8.8,  11.1,  8.3,  8.4,  9.6  degrees.  The  pattern 
here  consits  of  a  jump  in  angle  followed  by  a  smooth 
rise  to  the  next  jump  followed  by  a  smooth  rise.  This 
data  is  presented  as  1  part  of  Fig.  15.  The  dashed 
lines  mark  the  boundaries  of  the  indicated  cycle  number. 
The  constant  slope  is  the  RR  region  while  the  peak  is 
formed  by  the  RSR  ray.  From  Fig.  14  it  can  bo  seen 
that  the  peak  formed  by  the  13  cycle  RR  ray  does  not 
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continue  into  the  valley  formed  by  the  14  cycle  RSR 
ray  therefore,  there  is  a  large  chanqe  in  the  vertical 
angle  for  the  13  cycle  ray  as  shown  in  Fig.  15.  Prior 
to  750  km  this  peak  is  presented  in  each  ray  cycle. 

As  range  is  increased  beyond  1100  km  the  peak  formed 
by  the  N  cycle  RR  ray  reaches  further  into  the  valley 
of  the  N+l  cycle  RR  lay  and  ultimately  the  RR  rays 
dominate  the  vertical  angle  and  thus  the  bearing  bias. 

At  that  range  the  plot  similar  to  Fig.  15  should  consist 
of  only  the  small  slopes.  As  the  beam  reaches  broadside 
the  bias  in  bearing  estimation  decreases  rapidly. 

Unfortunately,  Figs.  14  and  15  are  not  as  readily 
constructed  for  the  cases  including  the  fronts.  The 
arrival  pattern  is  available  for  only  the  specific 
ranges  examined.  The  bias  in  bearing  estimate  is 
presented  for  1°  and  5°  beamwidths  in  Tables  1,  2  and 
3.  Fig.  If  presents  the  bearing  bias  at  endfire  for 
the  sampled  ranges.  From  750  km  to  850  km  a  pattern 
similar  to  that  in  Fig.  15  is  present  but  once  the 
fronts  are  encountered  the  data  is  not  sampled  finely 
enough  to  discern  a  new  pattern.  However,  it  can  be 
seen  that  the  first  front  does  decrease  the  bias  as 
predicted  by  the  focusing  effect. 

The  multipath  vertical  angle  pattern  for  the 
acoustics  passing  through  both  fronts  is  presented  in 
the  upper  portion  of  Fiq.  12.  It  can  be  seen  that  the 
number  of  arrivals  has  increased  and  there  are  several 
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rays  with  lower  vertical  angles  at  the  receiver  when 
compared  to  the  multipath  passing  through  the  first 
front  only.  Figure  16  includes  the  bearing  bias  at 
endfire  for  ranqes  into  the  2nd  front.  Although  the 
decrease  in  bias  is  not  as  smooth  as  at  the  first 
front  it  does  exist  indicating  a  continuing  focusing 
of  the  multipath  structure. 

We  will  now  briefly  examine  the  effects  of  the 
same  frontal  system  on  propagation  between  a  shallow 
source  in  warm  water  and  a  deep  receiver  in  cold  water. 

The  multipath  structure  as  a  function  of  range  is  shown 
in  Fig.  17.  A  repeating  pattern  is  present  up  to  1090 
km  indicating  the  passage  of  the  source  through  the 
first  front  (front  2  in  Figs.  1  and  16)  has  done  little 
to  affect  the  propagation.  However,  as  the  source  passes 
throuqh  the  second  front  (front  1  in  Figs.  1  and  16)  a 
siqnificant  change  in  the  multipath  structure  takes  place. 
The  major  arrival  pur  changes  its  vertical  angle  at  the 
receiver  from  3°  to  6°.  Thus  the  frontal  system  combined 
with  this  source  and  receivci  orientation  acts  as  a  weak 
dispersive  lens  once  the  source  is  in  the  warm  water. 
Figure  18  depicts  the  bias  in  the  bearing  estimate  at 
the  endfire  aspect  for  receivers  located  in  both  warm 
and  cold  water.  The  different  patterns  in  the  bias  for 
the  fo  cases  arc  present  from  ~b0  to  850  km.  Beyond 
the  second  front  the  bias  ha.;  decreased  5°  for  the  warm 
water  receiver  while  it  has  increased  3°  for  the  cold 
water  receiver. 
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V.  Results  and  Conclusions 

The  effect  of  an  ocean  front  on  the  multipath  vertical 
angles  and  the  associated  bias  in  bearing  estimation  was 
found  to  be  highly  dependent  of  array  orientation.  At  end- 
fire  the  maximum  bias  was  1 1(  while  beyond  30°  from  endfire 
the  bias  was  reduced  to  less  than  2°  for  the  various  ranges 
and  frontal  conditions  examined.  The  bias  was  found  to 
exhibit  a  somewhat  periodic  pattern  with  range  with  the 
amplitudes  of  the  variations  decreasing  with  increasing 
range.  Beamwidth  was  found  to  affect  the  bias  in  bearing 
estimate  to  a  lesser  degree  and  the  apparent  bearing  of  the 
source  was  always  directed  towards  broadside. 

It  was  also  found  that  the  frontal  system  acted  as  a 
convergent  or  dispersive  lens  depending  on  the  direction 
of  propagation.  For  a  warm  water  receiver  the  frontal  system 
decreased  the  bias  in  the  bearing  estimate  at  endfire  by 
5°  while  this  bias  was  increased  by  3°  for  a  cold  water 
receiver.  It  was  further  shown  that  a  narrow  beamwidth 
horizontal  array  should  be  able  to  separate  individual 
multipath  arrivals  at  endfire. 
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APPEND  T X 

Multipath  sound  propagation  through  a  cold  to  warm 

water  front,  subcripted  c  and  w  respectively,  will  undergo 

a  focusing  effect  when  the  ray  passes  under  the  anomaly. 

Referring  to  Fig.  Al,  the  ray  travels  from  cold  to 

warm  water,  passing  under  the  anomaly  of  depth  Zft.  Here 

Zv  and  Z  refer  to  the  vertex  depths  near  the  surface, 
w  c 
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